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ABSTRACT 

We investigate the change in ionizing photons in galaxies between 0.2 < 2 < 0.6 using the F2 
field of the SHELS complete galaxy redshift survey. We show, for the first time, that while the 
[O III]/H/3 and [O III]/[O II] ratios rise, the [N Il]/Ha and [S II]/Ha ratios fall significantly over the 
0.2 < z < 0.35 redshift range for stellar masses between 9.2 < log(M/M 0 ) < 10.6. The [O III]/H/3 
and [O III]/[O II] ratios continue to rise across the full 0.2 < z < 0.6 redshift range for stellar masses 
between 9.8 < log(M/M 0 ) < 10.0. We conclusively rule out AGN contamination, a changing ISM 
pressure, and a change in the hardness of the EUV radiation field as the cause of the change in 
the line ratios between 0.2 < z < 0.35. We find that the ionization parameter rises significantly 
with redshift (by 0.1 to 0.25 dex depending on the stellar mass of the sample). We show that the 
ionization parameter is strongly correlated with the fraction of young-to-old stars, as traced by the 
H/3 equivalent width. We discuss the implications of this result on higher redshift studies, and we 
consider the implications on the use of standard optical metallicity diagnostics at high redshift. 

Subject headings: galaxies: starburst—galaxies: abundances—galaxies: fundamental parameters 


1. INTRODUCTION 

An understanding of the conditions under which stars 
formed at all epochs is critical in constraining theo¬ 
retical models of galaxy evolution. Recent observa¬ 
tions of z > 1 star-forming galaxies suggest that the 
ISM state and/or the hardness of the EUV radiation 
field was more extreme in the past than in the present 
day. Galaxies at high redshift display larger line ratios 
([O III] /H/3 and/or [ N II] /Ha) than l o cal galaxies (e.g. , 
Hainline et al.ll2009t iBian et al.l 120101: lYa.be et all 120121 

Kewlev et al.1 l2013bl: iHolden et all 12014 ISteidel et al l 

2014 and references therein). These enlarged line ra¬ 
tios have been interpr eted in terms of a contribution 
from an AGN or shocks (iGroves et alll2006l: I Wright et al.1 
20101 [Trump et aD 120111 iForster Schreiber et al.l 12014 
Newman et al . 2 i 014), a larger nitrogen abundance 
(Mas ters et al.l 12014. a larger ionization parameter 
( Brinchmann et al.l 120081: IS teidel e t al.l 12014 1. and/or a 
higher ISM pressure dShirazi et al.1 l2014af) . Measure¬ 
ments of the ionization parameter and electron den¬ 
sity of the ISM using rest-frame optical line ratios sup¬ 
port the general picture of a larger ionization param¬ 
eter or electron density - altho u gh the useable sam¬ 
ples are small (iLiu et all 120081 Mainline et al.l 120091: 
Bian et al.1 120101 INakai ima et al.l l2012t iShirazi et al.1 
2014allbl : INakaiima fe Quchill20l41 . However, some high 
redshift gala xies have electron densities similar to local 
galaxies (see lRiebv et al.l[2011t iBavliss et all 1 2 01 4) . im- 
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p lying a different cause of the anomalous line ratios. 

iKewlev et all (j2013aT ) used theoretical simulations to 
show how the [N II] /Ha and [O III] /H/3 ratios of star 
forming galaxies and AGN may change with redshift, 
given four sets of extreme assumptions. If the ISM con¬ 
ditions change with redshift, then the standard optical 
classification schemes and other line ratio d iagnostics 
may n ot be applicable at high redshift (jKewlev et al.1 
MM)- Here, we investigate the change in optical line ra¬ 
tios across the intermediate redshift range 0.2 < z < 0.6 
using the F2 field of the SHELS galaxy redshift survey. 
The SHELS survey is 95% complete to a limiting mag¬ 
nitude of R = 20.6, allowing the optical line ratios to be 
investigated as a function of redshift and stellar mass. At 
high redshift, emission-line selection effects are difficult 
to remove. Galaxies with low [O III] or [N II] luminosities 
may be missing from high-redshift samples, potentially 
accounting for a large fraction_of the o bserve d change in 
line ratios with redshift ( Juneau et al. l; 2014 ). The high 
level of spectroscopic completeness of the SHELS sample 
avoids these issues and allows one to robustly constrain 
the change in line ratios with time, as well as the cause(s) 
of the change in line ratios. In this paper, we use the 
SHELS sample to show that the [O III] /H/3, [N II] /Ha, 
and [S II] /Ha ratios change systematically from z = 0.2 
to z = 0.6, and we demonstrate that this change is caused 
by a change in the number of ionizing photons per unit 
area within galaxies, as a function of redshift. Through¬ 
out this paper, we we use the flat A-dominated cosmology 
as measur ed by the 7 year WM AP experiment (h = 0.7, 
f l m = 0.3: lKomatsu et all 1201 111 . 

2. SAMPLE AND DERIVED QUANTITIES 

The SHELS survey is based on t wo of the Deep Lens 
Survey Fields (|Wittman et alJl2006D . We use the F2 fiel d 
from the SHELS survey, described in lGeller et al.l (120141) . 
The 4 square degree field contains 12,705 spectroscopic 
redshifts for galaxies with r<20.6. The survey is 95% 
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complete to this limit. 

The SHELS spectra were obtained by iGeller et al.1 
( 2005) usin g the Hectospec multi - fiber s pectrograph on 
the MMT (iFabricant, et al.l 119981 12005( 1 . The spectra 
cover the full optical wavelength range from 3700 - 9100 
A, at a spectral resolution of ~ 5A. This high spectral 
resolution is sufficient to resolve the [S II] doublet, al¬ 
lowing electron density estimates to be made by stacking 
in bins of stellar mas s and redshift . Ex ample individual 
spectra are shown in IGeller et al.l (12014H . 

We have selected a sub-sample of the SHELS catalog 
according to the following criteria: 

• A redshift in the range 0.2 < z < 0.6. The lower 
redshift limit avo ids aperture effects, which can be 
large for z < 0.2 dKewlev et alJl20Ml . 

• The 4000A break index, D„4000 < 1.5. This se¬ 
lection limits the sample to stellar ages < 1 Gyr, 
and is insensitiv e to the metalli city o f the stellar 
population (e.g.. iPoggianti fc Barbarolll997H . 

The D„4000 index is defined as the ratio of the flux 
(f„) in the 4000 - 4100A and 3850 - 3950A bands 
(jBalogh et al.lll999ll . The D n 4000 index has been shown 
to be less sensitive to reddening effects and contains 
less uncertainties from the _stellar c ontin uum than the 
original D4000 index (jBalogh et al.1 Il999h . The 4000A 
break is produced mainly by the heavy element absorp¬ 
tion in the atmospheres of the older stellar population- 
the hotter and younger stars produce a smooth contin¬ 
uum in this region of the spectrum. Thus, the D n 4000 
index provides a relative measure of the age of the stel¬ 
lar population. The D„4000 distrib ution for SHELS 
is des cri bed in IFabricant et al.l ((2008( 1 and IGeller et al.1 
1(2014(1 . IGeller et al.1 (12014H showed that the SHELS 
D„4000 distribution is bimodal, with D„4000 < 1.5 se¬ 
lecting for actively star-forming galaxies and D„4000 > 
1.5 selecting for quiescent galaxies. 

2.1. Stacking Method 

To avoid emission-line selection effects, we stack the 
data in bins of stellar mass, as described in IGeller et al.1 
(]2014f) . We divide our sample into bins of stellar mass 
in 0.2 dex increments. At least ~ 50 galaxies per bin are 
required to provide measurable [O III] and H/3 emission¬ 
lines in each stacked spectrum. The spectra and ob¬ 
servational uncertainties are linearly interpolated to a 
common rest-frame wavelength vector based on the red¬ 
shift of the bin. The interpolated rest-frame wavelength 
vector has a spectral resolution of 1.5A and spans an ob¬ 
served wavelength range of A3500 — 9100 A. The average 
flux at each resolution element is ta ken and the errors are 
added in quadrature. iGeller et al.l (1 2014 ) conducted an 
in-depth analysis of the use of stacking for the SHELS 
sample. They showed that stacking does not bias line 
ratio analysis for SHELS. 

2.2. Stellar Masses 

Stellar masses were derived for our sample from 
the broad-band photometry using the Le Phare code 
by Arnouts fc Ilbert ( privat e communication), with 
the iBruzual fc Charlotl ((20031 1 stellar templates as in¬ 
puts. These stellar templates were given for two 


metallicities and for seven exponentially decreasing 
star formation models (SFR oc e 4 / T ) with T = 
0.1, 0.3,1, 2,3, 5,10,15, 30 Gyr over a stellar population 
age between 0 — 1 3 Gyr. To cor r ect fo r dust attenua¬ 
tion, we apply the iGalzetti et ah ( 1200 0) law with E(B- 
V)= 0 — 0.6. The Le Phare code provides a mass prob¬ 
ability distribution. We use the median of the mass dis¬ 
tribution as the stellar mass of our galaxies. 


2.3. Ionization Parameter 

The ionization parameter, q , is defined as the number 
of hydrogen ionizing photons passing through a unit area 
per second divided by the number density of hydrogen 
atoms. The volume-weighted mean ionization parame¬ 
ter, q , can be defined in terms of the Stromgren radius 
R s for a filled spherical H II region as 


_ 2 2 ' 3 S t (t) 
^ 47 jR 2 s n 


( 1 ) 


where n is the total number density, and /?*(£) is the 
number of ionizing photons produced by the exciting 
stars. 

A larger ionization parameter at high redshift has been 
proposed as a pote ntial cause of the change in line ratios 
with redshift (e.g..lBrinchmann et al. 2008; Shiraz i et al.1 
I2014al ISteidel et al.ll2014H . However, most high redshift 
galaxy spectra do not contain sufficient emission-lines to 
reliably constrain the ionization parameter. Fortunately, 
we are able to constrain the ionization parameter for 
SHELS, and to investigate how the ionization parameter 
changes both as a function of redshift and as a function 
of stellar mass. 

We measure the ionization parameter using the 
[O III] / [O II] ratio using the theoret i cal ite rative ap¬ 
proach de scribed i n iKewlev fe Douital (|2002h . and up¬ 
dated in iKobulnickv fc Kewlevl ((20041) . We use the 
equivalent width ratio o f [O Illl/fO II] in t hese calcu¬ 
lation s. IKobulnickv fc Phillinsl (1200311 and iZahid et al.1 
(j2011h showed that the [O III] / [O II] flux ratios are well 
represented by the equivalent width ratio of [O III] / [O II]. 
For redshifts z < 0.4, we have verified that the use 
of equivalent widths yields the same results as using 
extinction-corrected line fluxes. 


3. AGN CONTAMINATION 

An AGN produces a hard ionizing radiation field which 
leads to elevated [N II] /Ha and [O III] /H/3 line ra¬ 
tios. A larger relative contribution from an AGN at 
high redshift has been proposed as t he cause of the rise 
in [O III] /H/3 ratio with redshift, (|Groves et al.1 120061 : 
iWright, et alJl20lfl : iTrumo et al.ll201ll ). Sample selection 
is likely to play a major role in whether AGN contribute 
to the observed line ratios in a given sample. 

The effect of residual AGN contamination on stacked 
data is unknown. However, the SHELS sample also al¬ 
lows us to investigate the effect of AGN on the [N II] /Ha 
and [O III]/H/3 optical line ratios ratios for galaxies be¬ 
tween 0.2 < z < 0.38 and on the [S II]/Ha ratio for 
galaxies between 0.2 < z < 0.35. We remove AGN and 
composites using the [O III] /H/3 versus [N II]/Ha op¬ 
tical di agnostic diagram, fo rming a purely star-forming 
sample (|Kewlev et al.l[200 6 j (Figure[l]top panel)). 
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In Figure [Q (middle and lower panels), we show how 
the optical line ratios change with redshift for each stel¬ 
lar mass bin (colored curves). The right panel shows the 
effect of AGN contamination on the line ratios in the 
stacked data. The hard ionizing radiation field from an 
AGN strongly affects the [O III]/H/3 ratio because more 
0+ ions are ionized into 0 ++ ions. It is clear that if AGN 
are not removed, the [O III]/H/3 ratio is strongly affected 
by AGN for stellar masses larger than logM/M Q > 10. 
Galaxies w ith l owe r stellar masses contain few AGN (e.g., 
iGroves et al.ll2006lf . Therefore, for z > 0.4 - where we 
cannot remove the AGN - we have restricted our analy¬ 
sis to galaxies with log(M/M Q ) < 10 to minimise AGN 
contamination. 

4. THE CHANGE IN LINE RATIOS WITH REDSHIFT 

Figure [2] shows that the [O III]/H/3 ratio rises sig¬ 
nificantly between 0.2 < z < 0.6. Although this ra¬ 
tio is also a strong function of stellar mass (Figure [2jr) 
nonetheless, the stellar mass does not appear to affect 
the size of the change in [O III] /H/3 shown in Fig¬ 
ure ca> for each mass bin. The [O III]/H/3 ratio ra¬ 
tio rises by 0.28 ± 0.01 dex across the 0.2 < z < 0.6 
redshift range for all stellar masses in the range 9.8 < 
log(M/M 0 ) < 10.0. This rise is consistent with high 
redshift studies which find large [O III]/H^ lin e ratios in 


star-forming galaxies with 2 : > 1.5 (Hainline et al. 

2009; 

Bian et al. 2010; Rigbv et al.i 20111 Yabe et al.l 

2012; 

Kewlev et al. 2013bt Holden et al. 20141 Steidel et al. 

2014f Sanders et al. 20lJ|). We will now briefly 

mves- 


tigate possible mechanisms which could produce this ob¬ 
served rise in the [O III]/H/3 ratio. 


4.1. Mean pressure in the HII regions 

Larger pressures in the H II regions of galaxies may 
cause a rise in the [O III]/H/3 ratio. This pressure can 
be inferred from the electron density of the gas. In a 
fully ionized plasma with an isobaric density distribu¬ 
tion, the particle density is defined in terms of the ratio 
of the mean ISM pressure, and mean electron temper¬ 
ature, T e , through n = and n is related to the 

electron density through n = 2n e (l + H e /H). For an 
ionized gas, the electron temperature is ~ 10 4 K and 
the electron densit y is t herefore propo rtional to the ISM 
pressure (see e.g.. iDopita et al.M2006l for a discussion). 
We can trace the electron density in the SHELS sam¬ 
ple for 2 < 0.35 using the [S II] A6717/[S II] A6731 ra¬ 
tio. We find that the difference in [S II] ratio between 
0.2 < z < 0.35 is consistent with no change with redshift; 
A([S II] A6717/[S II] A6731= 0.05 ± 0.05 from z = 0.2 to 
z = 0.35. 

Furthermore, a larger ISM pressure increases all three 
[N II], [S II], and [O III] ratios simultaneously because 
the larger pressure suppresses IR fine structure cooling, 
leading to more collisional excitations of all th ree [N II], 
[S II], and [O III] lines (see Figure 2 in iKewlev et al.1 
I2013af) . However, in the SHELS sample between 0.2 < 
z < 0.35, the log([N Il]/Ha) and log([S II]/Ha) ratios 
actually decrease with redshift by ~ 0.03—0.3 dex (panels 
c and g of Figure [1}. We therefore rule out increasing 
pressure with redshift as the cause of the observed change 
in SHELS line ratios across the 0.2 < z < 0.35 range. 
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Figure 1. (Upper panel) the optical diagnostic BPT diagram 
showing how AGN are removed from the SHELS 0.2 < 2 : < 0.38 
sample. Both AGN and composites are removed, yielding a purely 
star-forming sample. The middle panels show the the [N II]/Ha, 
and [O IIl]/H/31ine ratios versus redshift for the SHELS stacks be¬ 
tween 0.2 < 2 : < 0.38 where AGN have been removed (left panel) 
and AGN are included in the stacks (right panel). The lower panels 
give the [S II]/Ha ratio versus redshift for the SHELS stacks be¬ 
tween 0.2 < 2 : < 0.35 where AGN have been removed (left panel) 
and AGN are included in the stacks (right panel). The colored 
curves correspond to the stacked spectra in stellar mass ranges, as 
shown in panel d. The errors in the line ratios and the standard 
error of the redshift for each stack (shown) are smaller than the 
filled circles, thanks to the large number of galaxies in each stacked 
spectrum (> 500). AGN strongly affect the [O III]/H/5 ratios of 
the larger stellar mass stacks (log(M/MQ)> 10.0). 

4.2. Hardness of the radiation field 

Recently, a hard ionizing radiation field has been pro- 
posed as the cause of the large ]Q IIIl/H/3 line ratios in 
high redshift galaxies (jSteidel et al .1120 1 lL iStanwav et al.l 
I2014H . The effect of a hard ionizing radiation field on 
the optical line ratios is clearly illustrated in Figure [T] 
by comparing the panels with and without an AGN for 
high mass galaxies (log(M/Mg)> 10). A hard ionizing 
radiation field from either an AGN, shocks, or from a 










































4 



Figure 2. (a) The [O III]/H/3 ratio versus stellar mass for lines 

of constant redshift, as indicated in the legend, (b) the [0 III] /H/3 
ratio versus redshift for lines of constant stellar mass. The rela¬ 
tionship between colors and stellar mass is shown in panel (a). 
The [O III]/H/3 ratio rises by 0.2 — 0.3 dex between redshifts 
0.2 < z < 0.6. The change in [O III]/H/3 does not depend on 
stellar mass. 

hot stellar population (such as Wolf-Rayet stars) ionizes 
more neutral nitrogen and sulphur in to N + an d S + , as 
well as ionizing Q + into 0 ++ (e.g., IKewlev et all 120011 : 
iLevesoue et al.l[20Tof ). The [S II] emission-line is partic¬ 
ularly sensitive to the hardness of the ionizing radiation 
field. In a hard ionizing radiation field, a large, warm 
and partially ionized zone containing S + ion extends to 
the edge of the photoionised nebula, producing enhanced 
[S II] /Ha ratios. Because both [N II] /Ha and [S II]/Ha 
ratios fall with redshift, we can rule out a hard ionizing 
radiation as the cause of the change in SHELS line ratios 
with redshift. 

4.3. A changing ionization parameter with redshift 

A larger ionization parameter might produce the ob¬ 
served change in line ratios. In Figure [3] we show that 
the ionization-parameter sensitive line ratio [O III] / [O II] 
rises with redshift (by up to 0.25 dex). This change in 
[O III] / [O II] ratio corresponds to a change in the ion¬ 
ization parameter of up to Alog(g) ~ 0.25 cm s -1 . Al¬ 
though an increase in pres sure can lead to an apparent 
rise in q of this magnitude (iDooita et al.ll2014 ), the pres¬ 
sures needed would be in excess of log P/k > 6.OK cm -3 , 
which is excluded by an absence of any change in the elec¬ 
tron density as measured by the [S II] A6717/[S II] A6731 
ratio with redshift in the SHELS sample. 

Thus, any change in the hydrogen ionizing photon flux 
must be produced either by increasing the hardness (i.e. 
changing the slope) of the ionizing EUV radiation field, 
or by scaling up the whole ionizing radiation field by a 
some factor. We have already ruled out a change in the 
hardness of the ionizing radiation field. We now investi¬ 
gate the prospect of scaling the radiation field. 

A larger hydrogen ionizing photon flux on the inner 
edge of the nebula could either be produced by collective 
effects such as a larger number of stars within an individ¬ 
ual H II region, or by geometrical effects such as a gener¬ 
ally closer proximity of the stars to the ionisation front - 
as in blister H II regions. If the first hypothesis is correct, 
we would expect the ionization parameter to be related 
to the current star formation rate. This hypothesis can 



Figure 3. The ionization-parameter sensitive [O III] /[O II] ratio 
as a function of redshift for our stacked data in bins of stellar mass 
(colored curves). The [O III] /[O II] ratio rises with redshift for all 
stellar masses. 

be tested by comparing the [O IIl]/[0 II] ratio or the 
derived ionization parameter with the H/3 emission-line 
equivalent width, because the Balmer equivalent width 
traces the relative fract ion of young stars to old stars in a 
galaxy (see e.g.. iLeithereifl2005l for a review). Indeed, we 
find a strong correlation (Figure[4]) between these quanti¬ 
ties. The Spearman-rank correlation coefficient between 
the ionization parameter and the H/3 equivalent width 
is 0.81, and the probability of obtaining this value by 
chance is formally zero (6.6 x 10 -5 ). We conclude that a 
larger number of young to old stars within H II regions is 
likely to be principally responsible for the change in ion¬ 
ization parameter with redshift in the SHELS sample. 

If the ionization parameter increases with redshift, 
then this will have important implications on the use 
of line ratios to derive metallicity at high redshift. In 
particular, the [N II]/Ha ratio and the [O III] /H/3 ra¬ 
tios will be impacted strongly by a change in ionization 
parameter. Both of these line ratios are in widespread 
use in high-z metallicity studies due to their insensitiv¬ 
ity to extinction and flux calibration. If the ionization 
parameter is not taken into account, the [N II]/Ha ratio 
will underestimate the metallicity. The effect of ioniza¬ 
tion parameter on the [O III] /H/3 ratio is more compli¬ 
cated because the [O III] /H/3 ratio is double-valued with 
metallicity. At high metallicity (log(0/H) + 12> 8.4), 
the [O III] /H/3 metallicity will be underestimated, while 
at low metallicity (log(0/H) + 12< 8.4), the [O III]/H/3 
metallicity will be overestimated. We recommend that 
the ionization parameter and metallicity be derived si¬ 
multaneously to ensure reliable metallicity estimates for 
non-local galaxies. 

5. CONCLUSION 

We investigate the change in optical line ratios across 
0.2 < z < 0.6 for the SHELS survey. For 9.8 < 














5 



Figure 4. The [O III]/[O II]— ratio (upper panel) and the ioniza¬ 
tion parameter (log q; lower panel) versus the H (3 equivalent width 
for our SHELS sample in bins of stellar mass (colors as in Figure 1) 
and in bins of redshift (symbols) 

log(M/M 0 ) < 10.0, the [O III]/H/3 and [O III]/[0 II] 
ratios rise (by ~ 0.2 dex) across the full 0.2 < z < 0.6 
range. For 0.2 < z < 0.35 we show that while the 
[O III]/H/3 and [O III]/[0 II] ratios rise, the [S II]/Ha 
and [N II]/Ha ratios fall significantly (by 0.05 — 0.3 dex) 
for stellar masses betweeon 9.2 < log(M/M 0 ) < 10.6. 
We explore these changes in terms of the electron den¬ 
sity (or ISM pressure) of the gas, the slope of the ionizing 
radiation field, and the relative fraction of young to old 
stars. We rule out a change in the electron density and 
a change in the hardness of the ionizing radiation field 
as dominant causes of our observed change in line ratios 
with redshift. 

We examine the effect of AGN contamination. 
We show that AGN substantially contaminates the 
[O III]/H/3 line ratios at stellar masses (log(M/M 0 ) > 
10). AGN removal is critical for high-z studies that inter¬ 
pret the [O III] /H/3 ratio in terms of star-forming galaxy 
properties. However, contribution by AGN does not af¬ 
fect the general trends observed in the SHELS [N II]/Ha 
and [S II] /Ha line ratios. It is unclear whether this lack 
of AGN sensitivity will hold at high-z, given potentially 
different ISM conditions and AGN properties. 


We show that the observed change in SHELS line ratios 
across 0.2 < z < 0.35 is dominated by a rise in the ioniza¬ 
tion parameter with redshift. The geometry, stellar mass 
contained, and the age of the stellar populations within 
H II regions may all contribute to the observed ionization 
parameter. The observed change in ionization parameter 
is strongly correlated with the fraction of young to old 
stars, as traced by the H/3 equivalent width. 

We emphasize that while the ionization parameter 
dominates the change in line ratios observed in SHELS, 
the ionization parameter may not be the dominant cause 
of a change in line ratios for different redshift ranges or 
different stellar mass ranges, particularly where shocks 
or AGN contamination may be present. 
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